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ABSTRACT 


Several earlier investigators who have considered hypersonic gas dynamic 
heating have suggested the use of empirical correlation formulas such as: 



where Q is the stagnation point heat transfer rate, R is the vehicle nose radius, 
w 

n is the ambient density and is the flight speed of the vehicle. 

' 00 

However, since the magnitude of the coefficient C depends on the contribu- 
tion of the actual chemical species present in a given planetary atmosphere, C is 
not known "a priori" but can be obtained theoretically by solving the multicomponent 
boundary layer equations for a variety of foreign gases, which has not been accom- 
plished heretofore. 

The purpose of the present investigation, therefore, was to determine the 
thermochemical effects of foreign planetary atmospheres upon hypersonic stagna- 
tion region laminar heat transfer and skin friction, and to attempt to obtain, if possi- 
ble, a simple correlation of the results with molecular weight, over a range of free 
stream molecular weight from 2 to 44. 
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INTRODUCTION 


With the advent of space exploration, interest is growing in the problem of 
entry into planetary atmospheres other than Earth. Thus, the problem of predicting 
the aerodynamic heat transfer rate becomes complicated by the presence of chemi- 
cal species other than those present in high temperature air. 

The theoretical investigations of the determination of heating rates in disso- 
ciating air flows are well documented and some representative studies concerning this 
problem area appear in References 1 to 5. A comparison (Ref. 5) between these 
various theoretical predictions of stagnation point heating in dissociated air indicates 
good agreement with a maximum deviation of approximately ± 20 per cent from a 
mean value. 

The problem of predicting the stagnation point heat transfer in gases other 
than air has led some of the earlier (Ref. 6-10) investigators to suggest the use of 
generalized correlation formulas such as: 



In the above equation, the value of the coefficient C is intended to reflect the contri- 
butions of the various chemical species present in any atmosphere to the transport 
and thermodynamic properties of the resulting dissociated gas mixture. 

Gazley (Ref. 6) estimates the coefficient C to be proportional to the square 
root of the ratio of the free stream Mach number to the Reynolds’ number for the 
laminar continuum regime of planetary entry. This calculation leads to values for C 
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which are similar for Earth, Mars and Venus to within ± 3 per cent. Chapman 
(Ref. 7) relies on the dissociated air theory of Lees (Ref. 2) to obtain estimates of 
heating rates in gases other than air and finds that C ~ Pr^"^ jli (y - l/y)^^. 

In another study, Lees (Ref. 8) offers an equation for engineering purposes in which 
the constant C is to be taken the same for the planets Earth, Mars and Venus. In 
his correlation equation, however, provision is made for the variation of surface 
temperature. 

The first rigorous theoretical investigation of hypersonic heat transfer in 
a gas different from air was that performed by Hoshizaki (Ref. 10). In his study, 
Hoshizaki calculates the equilibrium boundary layer flow at the stagnation point of a 
blunt body in carbon dioxide gas as well as in air. Preferential diffusion of the 
various chemical species is neglected, and the effects of dissociation are taken into 
account by means of the "total" thermodynamic and transport property concept, 

(Ref. 11, 12, 13). Hoshizaki's correlated heating rates in air compare favorably 
with the earlier theoretical studies of Lees (Ref. 8) and Chapman (Ref. 7). Further- 
more, he finds that the coefficient C for CO is approximately 10 per cent higher 
than that for air. 

Experimental efforts to determine heat transfer in gases other than air 
include hypersonic shock tube (Ref. 14, 15) and ballistic range (Ref. 16) measure- 
ments in pure nitrogen, pure carbon dioxide and mixtures of nitrogen, carbon dioxide 
and argon. 

The intent of this present study is to consider the thermochemical effects 
rigorously and to determine theoretically the aerodynamic heat transfer rate and skin 
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friction at the surface of a vehicle traveling at hypersonic speed in a planetary atmos- 
phere of arbitrary chemical composition. In order to accomplish this end, the laminar 
heating rates and surface shear stresses in eight dissociating gas mixtures having 
various degrees of chemical complexity are first determined. Then, simple correla- 
tion equations based on the molecular weights of the undissociated gases are derived. 

The gaseous systems considered in this study are as follows: three diatomic 
gases including N 2 , and C> 2 ; a triatomic gas, CC> 2 ; and several mixtures including 
four component "Air", i.e. , O, N and N g (simulating Earth's atmosphere), and 
several six component mixtures of the dissociation products of N and CO , i.e. , 

C0 2 , CO, 0 2 , O, N 2 and N (simulating the Martian and Venusian atmospheres). 

The approach used in this study consists of first selecting a suitable gas model 
based on the predominant species existing in the particular planetary atmosphere under 
investigation. The next step is the determination of the thermochemical properties 
of the gas existing in the flow field between the detached shock wave and the entry 
vehicle, by solving the normal shock and isentropic stagnation relations assuming ther- 
mochemical equilibrium. Finally, the theoretical laminar heat transfer rates are 
determined from the exact solutions of the multicomponent laminar boundary layer 
equations, as is the skin friction coefficient. 

In this study it is necessary to solve the normal shock equations and isentropic 
stagnation equations in order to provide self-consistent outer boundary conditions for 
the boundary layer computations. Also, by this means, it is possible to relate une- 
quivocally the boundary layer solutions to the external ambient conditions for each of 
the gas systems under consideration. 
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The appropriate governing boundary layer equations for a compressible viscous 
multicomponent chemically reacting gas are presented, and the transport and thermo- 
dynamic properties of the gaseous mixture are evaluated as a function of the equilibrium 
gas composition and temperature. In the numerical integration of the boundary layer 
equations, the PrandU number of the mixture and the multicomponent Lewis numbers 
are evaluated at the surface conditions. The viscosity of the mixture and all thermo- 
dynamic coefficients are evaluated locally throughout the viscous layer. The surface 
and outer edge boundary conditons are satisfied, assuming that the gas is m a state 
of local thermochemical equilibrium throughout the boundary layer. Numerical solu- 
tions are obtained on an IBM 7090 digital computer, and the results presented m the 
form of correlation equations and graphs. The calculations cover an appropriate range 
of flight speeds (i. e. , approximately 10,000 to 50,000 ft./sec. ) necessary to correlate 

the dissociation effects for each of the gaseous systems considered. 

In addition to the correlated heat transfer rate, the skin friction coefficients 
are also correlated, and are presented for the gaseous systems treated. Comparisons 
are made between the shear stress and heat transfer, and the agreement with Reynolds' 
analogy is noted, Finally, the new theoretical results are compared with other avail- 
able theoretical predictions and available experimental data. The agreement with 
other theoretical results and with experimental data is found to be good. 
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thermochemistry 


JCHtlVli&iiVi 

TI^Tconsider the determination of the thermochem.cal state of the gas 

surrounding a space vehicle during planetary entry. During hypersonic flight, m a 

. .. . (ree stream gas is dissociated and perhaps partially ionised 

planetary atmosphere, the free stream gas 

by the shock wave which precedes the vehicle. For known ambient con . ions in 
tree stream gas (f.e. pressure, temperature, density and chemical composition, an 
tor a given flight speed, one is able to determine the downstream conditions behin 
shock wave by a simultaneous solution of the Rankine-Hugoniot relations together wa 

an appropriate equation of state. 

These algebraic equations may be written: 


P„ v * = P 2 V 2 

P* + Poo Vjf = P 2 + P 2 V 2 2 


V' 


V r 


( 2 ) 

(3) 

(4) 


K + — = h 2 + 2 

tinn of mass momentum and energy across a normal shock 
and represent the conservation of mass, man 

wave, respectively. Equations 2, 3, and 4 relate the density, the gas v 
pressure and the enthalpy at the upstream location .to the downstream location 2 
behind the shock. In addition to the above conservation equations, one must intio uce 
M appropriate equation of state for the gas mixture in order to solve this system of 

equations. 

^ rp 

p = p — T 

The chief difficulty in solving this set of equations resides in the fact that in 
a chemically reacting gas, the enthalpy and the mean molecular weight are functions 
„f the gas composition. The composition of-this gas is in turn a function of the un own 
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temperature and pressure, and hence an iterative solution is required. Therefore f„ 
order to proceed, one requires a taiowledge of the equilibria of the resulting high 
temperature gaseous system downstream of the shock wave. 

In general, the coupling of the high temperature equilibrium thermodynamic 

properties of the gaseous system under study to the normal shock equations may be 

done in t wo ways. The most popuiar numerical procedure ffor norma, shock waves in 

air, Ref. 17-20, has been to store tabular data in die form of Mollier charts for the 

ranges of temperature, pressure, density, enthalpy, etc. iikely to be encountered 

The alternate procedure is to compute the necessary thermodynamic properties of the 

dissociated gas mixture as part of the numerical solution of the overall problem, 

(Ref. 21 and 22). 

The second of these methods is the procedure which was used in this study. 
This method, while being somewhat restrictive in terms of the number of chemical 
species one may conveniently handle numerically, is extremely flexible in that the 
relative abundance of any given chemical element in a particular system may be varied 
parametrically by simply changing one or more constants. 

Consider a planetary abnosphere in which the ambient undisturbed gas ahead 
of the shock wave is characterized, for example, by a mixture of carbon dioxide 
molecuiar nitrogen and argon. The number and type of chemical species present* 
the flow fieid between the shock wave and the entering space vehicle will, of course, 
depend upon the shock strength. For flight velocities sufficiently high to cause 
dissociation of the gas, one may expect the foliowing chemica. species to be present: 

C0 2 , CO, C, C 2 , O, 0 2 , N, N 2 , A, NO, CN, etc. While it is not practical, nor 
indeed necessary, to consider all of the possible chemical species (since some appear 
in oniy trace amounts,, one can discuss, ,» general, an „ component system and the 
n independent relationships necessary for the solution of the mathematical system. 
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In order to determine the composition of an n component chemically reacting 
gas mixture originating from k chemical elements, one has available k-1 relationships 
for the conservation of chemical elements, one relationship for the overall conservation 
of species, and n-k relationships dealing with the rate at which the chemical reactions 
are progressing. We begin with relationships for the conservation of chemical elements, 
which, for our purposes, take the form of mass balance ratios. 

For this system composed of four chemical elements (i. e. oxygen, nitrogen, 
carbon and argon) one may write three mass balance ratios as follows: 



( 6 ) 


In this notation, the mass of the chemical element (oxygen for example) is defined in 
terms of the mass fractions in the following manner. 


O = C 


O 


"Or 


+ ^o L 


i \7??i/0 °2j^iWj/0 2 


(7) 


In equation 7 the first two mass fractions account for the atomic and molecular oxygen 
in the free state. The first summation quantity includes all the atomic oxygen in a 
chemically bound (i. e. CO, NO) state, and the second summation quantity includes all 
of the chemically bound molecular oxygen. The masses of the other chemical elements 
are defined in a similar manner. 

The values of the three constants k^, k 2 , and k 3 are given by the relative 
abundance of the chemical species in the undissociated ambient atmosphere. As an 


example, k 2 is evaluated below. 

~ (>?o />Tco ) c co 2 „ 

k,=-g« % — (8) 

2 N °N 2 „ 

The relationship for the overall conservation of species takes the form of 
Dalton's Law for a static system in thermochemical equilibrium, and provides us 
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with the condition that the sum of all the mole fractions must be unity. 


L X = 1.0 (9) 

i 1 

th 

where i is the i chemical species in the n component mixture. 

For the case where the gas has equilibrated itself both thermally and chemi- 
cally, the remaining n-k relationships are available in the form of the equilibrium 
functions. The assumption of thermochemical equilibrium implies that the character- 
istic chemical reaction times are several orders of magnitude smaller than the available 
flow times and that the dissociation and recombination rates of the various chemical 
reactions are effectively infinitely fast. An example is shown below which considers 
the chemical reaction involving carbon monoxide, atomic oxygen and atomic gaseous 
carbon. 


C(g) + CO 

x co 

K Pco = X C(g) x o p 


(10) 


In the above, P is the pressure of the gas mixture, expressed in atmospheres. The 
Kp functions themselves are exponentially dependent upon the temperature, and for 
large ranges in temperature may be expressed in a simple Arrhenius form. The 
auxiliary relationships for the system are the expressions for the mass fractions and 
mean molecular weight: 


X.'ZJ? 

c. = = £ x.^ (11) 

1 "75? i 1 1 

After the composition of the gas has been determined at a given pressure and temper- 
ature, the enthalpy of the gas mixture may be calculated by means of the relationship: 

h = I C h (12) 

i 1 1 
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The enthalpy of a pure species, h., is given as: 



o 


dT + Ah 


fi 


(13) 


where A hr is the enthalpy of formation at the reference temperature, 
n 


It is noted that the resulting system of n algebraic equations is highly coupled 
and extremely non-linear and, in general, for n > 2, iterative techniques must be 
employed in order to obtain a solution to the system. 

We will now apply these general equations to the specific gas systems under 
consideration. The pure diatomic gas systems (i.e. oxygen, nitrogen, and hydrogen) 
each contain a single chemical element and hence the mass balance ratios (equations 6) 
are not required. Furthermore, the atomic and molecular species in equilibrium are 
directly determined by utilizing their K p relationship. Therefore, the composition of 
each diatomic gas system is determined by equation 9 together with the appropriate 
Kp function. The equilibrium reactions for the various diatomic systems are as 
follows: 

Diatomic oxygen system: 


0 2? =± 20 



x o p 


X 


o r 


(14) 


Diatomic nitrogen system: 


N 


2 *~ 


— 2N 


K. 


X N P 


N 


X 


N r 


(15) 


Diatomic hydrogen system: 



2H 


X H P 


X 


H, 


(16) 
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The triatomic gas (CO ) considered here contains two chemical elements and 

Li 

requires a singe mass balance ratio as follows: 

(^o 2 /Vco 2 ) C CO, 


o 


32 

12 


(17) 


1 = c = (^c/**co 2 ) c co 2 

In the dissociated mixture, we will consider a model consisting of the five species, 

CO , CO, O , O and C (g). Therefore, in addition to equations 9 and 17, three 
2 * 9 2 

equilibrium relationships are needed to complete the system. Two of these are given 


by equations 10 and 14 and the third is shown below: 
CO + O — C0 o 


X 


CO r 


( 18 ) 


K, 


co 2 


X CO x o P 


In the analysis of heat transfer phenomena, high temperature dissociated air 
is closely approximated by a four component mixture of atomic and molecular oxygen 
and nitrogen. The single mass balance ratio required for this system is: 


„ _ 0_ C °* C °2 
2 N C N + C N 0 


= 0. 3068 


(19) 


Equation 19 states that the composition of equilibrium air is 23.48 % oxygen and 76.529c 
nitrogen by mass in the free stream. Dalton's Law (equation 9), is, of course, used 
for all the gas systems. The two remaining equilibrium functions for the "Air" system 


are given by equations 14 and 15. 

Six chemical species are treated in the system for the C-O-N gaseous mixture, 

including CO , CO, O., O, N and N. Of the two mass balance ratios required for this 
2 2 ^ 


system, one is given by equation 17 and the other is shown below. 

♦fo/*co) c co + K>co 2 ) c co f 


c o + c o. 


k 2 = 


C N + C N, 


( 20 ) 
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In the C-O-N system, the numerical value of k 2 was altered to simulate three mixtures, 
corresponding to 90% COg - 10% N^, 50% C0 2 - 50% Ng and 10% COg - 90% Ng. The 
equilibrium relationships for this system are given by equations 14, 15 and 18. 

Having chosen a particular system for study together with a set of ambient 
free-stream conditions, the state of the gas is first determined at point 2 behind the 
normal shock by solution of equations 2 through 5. Then, an isentropic calculation 
(see reference 22) is performed along the stagnation stream line. A compatible set of 
boundary conditions at the outer edge of the boundary layer is thus determined. 


BOUNDARY LAYER CONSIDERATIONS 

The non-linear partial differential equations of change for a multicomponent 
chemically reacting gas are derived, for example, in Ref. 23 and include the conser- 
vation of mass, chemical species, momentum and energy as shown below: 


+ V • (py ) = O 

dp i 

— * v ■ < p i Xi> = *i 

dv 

p tt = V • 7T + L p. F. 
d t ~ i i — * i 

p 4 r = - V* Q + 77 :Vv + £ p. V. • F. 
dt - i 1-1 - 1 


( 21 ) 

( 22 ) 

(23) 

(24) 
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In the absence of thermal diffusion, the conservation equation for species i becomes: 


pm, 


where 


/ 3B. BB.\ 

: i( u 'dT + v bj 

Bj - Xj/% 


= IS 

,j*i -?n 2 ij s y . 


= w. 


(26) 


The stream wise component of momentum becomes: 


_ 3u _ Bu 
pu — + Pv — 


ap . a ... aii . 

3x By ^ 3 v ^ 


The normal component of momentum is: 
BP 


Sy 


= o 


The energy equation becomes: 

— BT B T 

PC ( u - — + v t — ) = u 
p ' Bx By' 


fi + M<&> + ( K l?) 


' By By 

E C .2. pa. -S- S w.h. 

i P,- ^y ^2 ij By i l l 


Upon introducing the Mangler-Dorodnitsyn (Ref. 2) transformations 


T) = 


P e U e 

yfH 


r y 6 p , 

/ r d y 

«' o P Q 


t 2 6, 

C = f puurdx 
J w 'w e 


(27) 


(28) 


(29) 


(30) 

(31) 


and assuming that local similarity holds, then equations 25 through 29 may be reduced 
to a set of ordinary non-linear differential equations. The diffusion equation for the i^ 
chemical species becomes: 


(" — L £ 
[ Pr j ti 




M* 2 


J- L..X. 
1J J, 


18 


w. 

1 




* m 


(32) 


where ^ = 2 dlnu g /dln £ and the stagnation point velocity gradient is assumed to be 
Newtonian. 


The conservation of momentum becomes: 

P 

+ f f + 0 

-v-i v i 


(' tf ¥?) 

\ / v 


V r l 


e ^2 


= 0 


(33) 
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The energy equation becomes: 



The thermodynamic properties and transport coefficients required for defini- 
tion of the physical problem include a specific heat, an enthalpy, a coefficient of 


viscosity and a self diffusion coefficient for each of the n pure species present in the 

2 _ 

boundary layer. In addition, — — 2- symmetric binary diffusion coefficents are 


required to represent the transport interactions between unlike pairs of particles. 


The transport properties may be calculated from the following equations 


(Ref. 23) 


1 /2 

(ff M.-ft T) ' 


“i = 

16 77a 2 S2.< 2 ’ 2 )* 
1 1 

(35) 

£) 

(2Tr^ 3 T 3 /M..) 1,/2 
3 i y 

(36) 

ij 

" 16 p hoJqJ 1 - 1 )* 



which are the rigorous kinetic theory formulae for the viscosity and the diffusion 


coefficients respectively. In order to evaluate these properties one requires a 

knowledge of the collision diameter a, and the reduced collision integral Cl*. 

In the above equations, M.. is the reduced mass of the binary interaction and a. . is the 

^ J 1 

simple average of the collision diameters for the pure species i and j. A more 
complete discussion of the transport coefficients is given in the appendix. 

With the exception of three chemical species (i. e. C(g), H, H^) the numerical 
values for the transport and thermodynamic properties used in this study are given in 
references 25 and 26. The thermodynamic properties for carbon gas were obtained 
from reference 27, while reference 28 supplied the thermodynamic data for the 
hydrogen system. The transport coefficents for these three species were calculated 
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using equations 35 and 36. The interactions involving carbon gas were assumed to 
follow a rigid sphere potential and the collision diameter for this atom was taken to be 
3.42 A. The Leonard-Jones 6-12 potential function was assumed applicable for the 
hydrogen system. The potential parameter (i. e. c/k) for this binary encounter 
(i.e. H-H„) is 38°K and the collision diameters of the atomic and molecular species 
are 3.38 A respectively. (Ref. 23, 29). To summarize then, with the exception of 
hydrogen, all molecule -molecule interactions are based on the use of the Lennard- 
Jones 6:12 potential, whereas all atom-atom and atom-molecule interactions are based 
on a rigid sphere potential. 

The boundary conditions on the velocity and temperature as applied to the non- 
dimensional boundary layer equations (eq. . 32-34) are now written. At the surface, 
the "no slip" condition gives a zero tangential component of velocity. 

u 


In the absence of mass injection or suction, the normal component of velocity 
vanishes, hence, • 


(37) 


-f 

w r 


w 


du 


1 l/2 


= 0 


P M ~a~ / & 
w w dx 


and the non-dimensional surface temperature becomes 


(38) 


The boundary conditions at the outer edge of the boundary layer for velocity and 
temperature are given by: 


(39) 


lim f = lim 0 = 1.0 

V (40) 

ri - °° T) -> 00 
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In addition to these usual boundary conditions imposed on the temperature and 
velocity, one must be concerned with the proper chemical constraints on this system of 
equations. In the ensuing discussion, it will be assumed that the local state of the gas 
composition within the viscous layer is one of thermochemical equilibrium. The 
analysis follows the same general format as given in the preceding section. For genei 

ality, let us consider an n component C-O-N gas system. 

At the outer edge of the boundary layer the temperature, pressure, enthalpy, 
etc. are determined by the methods described earlier (i.e. a normal shock computation 
followed by an isentropic compression). The chemical composition at the edge of the 
viscous layer is, of course, available as a result of this solution. 

Now, within the viscous layer, where diffusion effects cannot be neglected, the 
lighter-weight particles tend to diffuse faster than the heavier ones which results in 
preferential diffusion. Because of this, the k-1 relationships for the conservation of 
chemical elements will not be in the form of mass balance ratios (eq. 6). Instead, 
constraints are placed on the chemical source terms (w^ appearing in equations 32 and 
34. Since the masses of the indivudual chemical elements must be conserved, one may 

write for the C-O-N system : 

W = O; W = O (41) 

O N 

The conservation of mass of the element oxygen is written as follows: 



The interpretation of the above expression is similar to that of equation 7. 

For the overall conservation of mass, one has 

L W. - O (43) 

i i 

It is obvious that equations 41 and 43 together imply that the mass rate of formation 
of the element carbon is also zero. 
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To complete the analysis of the gas phase chemistry, it is noted that n - k re- 
lationships are still needed which specify the rates of the chemical reactions. From the 
assumption of local thermochemical equilibrium, one has the necessary information from 
the n - k equilibrium functions, in the form of equation 10. 

In this study, it was not found necessary to consider either catalytic or non- 
catalytic effects of the surface on the boundary layer chemistry because of the initial 
assumption of local thermochemical equilibrium through the gas layer. Although real 
engineering surfaces have finite catalytic efficiency, since the gas is already assumed 
to be in an equilibrium state, the chemical activity of the surface can in no way influence 
the gas state without producing a departure from equilibrium in the gas phase. It is 
usually under conditions of low density (i. e. , the low Reynolds number regime) where 
the influence of finite catalytic efficiency is observed, because in this case, the free 
stream gas is not initially in a state of thermochemical equilibrium as it enters the low 
density viscous layer. Only the high Reynolds number problem has been treated here. 

In the absence of mass transfer (transpiration cooling or ablation), the net 
mass flux at the surface must be zero. This means that the net mass flux of the k 
chemical elements must separately vanish at the surface. One is free to soecify k-1 
such relationships as follows: 



The net mass flux of the element oxygen for example, is defined to be: 

( 45 ) 

To complete the system, one makes use of Dalton's Law (eq. 9) and the n - k equi- 
librium K functions. 

P 
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DISCUSSIONS OF NUMERICAL RESULTS 


a) Thermochemistry and Boundary Layer Profiles. 

Comparisons of some of the thermochemical properties and compositions of 
the various gas systems are presented in figure 2 through 5 and in Table I. Figure 2 
shows the enthalpies of the pure cnemical species as a function of temperature. The 
reference temperature of 536. 7°R has been selected for these calculations, and the 
enthalpies of H , N , and O molecules in their standard state are taken to be zero at 

Li £ j 

this temperature, in agreement with customary thermodynamical definitions. The 
enthalpy of the hydrogen atom (not included in figure 2) was based on a heat of 
formation of 93,762 BTU/LB. and a specific heat of 4.97 BTU/LB. °R. It is also 
noted that the enthalpies of formation of carbon monoxide and carbon dioxide are 
negative with respect to the major compo.ients of air, (see Table II). 

The equilibrium dissociation levels of the three diatomic gases treated are 
shown in figure 3 as a function of temperature at a pressure level equal to one 
standard earth atmosphere. These were based on the equilibrium constants given in 
Table III. A comparison of the gases indicates that the equilibrium dissociation levels 
of oxygen and hydrogen are nearly the same function of temperature. Nitrogen, on the 
other hand, has a dissociation energy per mole approximately twice that of oxygen and 
hydrogen and consequently its dissociation occurs at correspondingly high gas tem- 
peratures. 

The equilibrium chemistry in terms of mole fractions of the carbon dioxide 
gas system is seen in figure 4. Three regimes of chemical activity are apparent. At 
the lower temperatures, carbon dioxide is dissociating to form carbon monoxide and 
molecular oxygen in accordance with the following reaction: 

CC> 2 ^=^ CO + 1/2 0 2 
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At slightly higher temperatures, and before the above reaction has gone to completion, 
the oxygen molecule dissociates, (i.e. 0 2 -*20). At the highest temperatures, the 
carbon monoxide breaks down to form atomic gaseous carbon and additional atomic 
oxygen. The mean molecular weight of this system therefore ranges from its low 
temperature value of 44 LB/ LB. -mole to 14. 67 LB/LB-mole for complete dissociation. 
The magnitude of the enthalpy of this gas ranges from -3850 BTU/LB at the reference 
temperature to over 17, 000 BTU/LB at a temperature of 16,000 °R. The resulting 
overall enthalpy change is nearly 21, 000 BTU/LB which is approximately twice the 
enthalpy change for air over the same range of temperature. This would imply that in 
a carbon dioxide atmosphere, higher flight velocities may be sustained (than in air, 
say) before appreciable ionization levels would be encountered. 

Figure 5 shows two gas models which may be used to simulate the air system. 
The four component model includes molecular and atomic oxygen and nitrogen, while 
the five component model introduces the additonal species nitric oxide. It is found 
that the four component model shown in figure 5 is entirely adequate for simulating 
the energy content of high temperature air over a wide range of altitudes and flight 
speeds. 

The normal shock and isentropic compression calculations provide another 
means for comparison of the different gas systems. Numerical results of a represent- 
ative set of normal shock and stagnation point computations are shown in Table I. 

Note that the three ambient conditions (i.e. flight velocity, temperature, and pressure) 
are fixed for all of the gas systems. The free stream density, therefore, increases 
as the mean molecular weight of the gas system. The free stream enthalpy is seen to 
be a direct function of the enthalpies of formation of the gases comprising each system. 

Consider next in Table I the temperature (T 2 ) immediately behind the normal 
shock. Since the hydrogen system has a large heat capacity and an enormous enthalpy 
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of formation for atomic hydrogen, it is not surprising to find a low shock temperature 
with an associated low dissociation level. The nitrogen, air and carbon dioxide systems 
are found to have very nearly equivalent shock temperatures. However, the oxygen 
system, due to the relatively lower dissociation energy of the oxygen molecule, is 
found to be completely dissociated for this flight condition, and consequently gives rise 
to a higher shock temperature. The shock pressure P^ is found to increase with 
increasing molecular weight of the particular gas system. This is expected from the 
free stream density variation with molecular weight. Referring to equation 3 and 
utilizing the hypersonic approximation, one has 


P V 


and then introducing the equation of state (5), one obtains 


p Ko 

P ~ — QQ 

2 UC T 


V 

00 


2 


(46) 


(47) 


Hence, for the same initial free stream conditions of pressure P , and temperature 
T^, and flight speed V m , the shock pressure P ^ is directly proportional to the mean 
molecular weight of the ambient gas. For those gases for which T 0 is nearly the 
same, the density p behind the normal shock has, in general, the same variation with 
free stream molecular weight as does the shock pressure. However, we note that since 
the oxygen system is completely dissociated, its gas density is slightly lower than that 
of the air or nitrogen systems. One notes further that the enthalpies h ^ of hydrogen, 
nitrogen, air and oxygen are approximately the same. However, the enthalpy of the 
carbon dioxide system is substantially lower due to its negative heat of formation. The 
absolute level of the enthalpies (i. e. h , h 0 , h ) for all gas systems will be dependent 
upon the selection of a particular reference temperature. However, the enthalpy changes 
and the other thermodynamic quantities (i. e. T, P. p,7)t.) are quite independent of any 
arbitrary selection of the reference temperature. 
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The stagnation region values of temperature, pressure, density and enthalpy 
are seen to increase slightly over their normal shock values. However, the com- 
position of the gas, as reflected by its mean molecular weight, has not changed 
significantly during the isentropic compression process. The stagnation point 
velocity gradient, given in Table I, was calculated following the modified Newtonian 
theory. The percent dissociation tabulated for each of the gases was based on the 
minimum mean molecular weight possible for each dissociated system and was 
calculated by the following equation: 

% Dissociation = — x 100 (48) 

If one had chosen to compare the various gas systems on the basis of a 
constant free stream density instead of a constant free stream pressure as indicated 
in Table I, then equation 46 shows that the stagnation pressures of all the gases will be 
identical for the flight speed considered. Equation 4 indicates that the stagnation 
enthalpies would be the same as those listed in Table I. Examination of equation 12 
and 13 and figure 2 reveals that the equilibrium shock temperature of the carbon 
dioxide mixtures will be less than that for air, oxygen or nitrogen. However, ex- 
amination of the equation of state for these conditions reveals that the shock density 
of pure carbon dioxide will still be greater than that for air. 

In general, the state of a gas stream is completely defined by two independent 
thermodynamic state quantities in addition to the macroscopic stream velocity. However, 
as will be shown later, one may predict the approximate aerodynamic heating rate in 
any arbitrary gas mixture by specifying only p^, and M. . 

The information contained in figure 6 Drovides a means of estimating an 
upper limit on the validity of the gas models considered in this study. The curves of 
free stream densities versus flight speed indicate a bound on the flight corridor in 
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which the gases shown are approximately less than one percent ionized, or, to the same 
order of approximation, less than ninety-nine percent dissociated. Here, we have 
definite evidence that a carbon dioxide atmosphere can sustain relatively higher flight 
speeds than the other gases shown before encountering appreciable ionization phenomena. 
In a pure hydrogen atmosphere, of course, the vehicle can travel at flight velocities in 
the range of 60 - 70 thousand feet per second before producing dissociation inside the 
shock layer. It should be noted that the free stream density and flight velocity of one 
of the cases that has been treated (i. e. oxygen at 24, 000 ft. /sec. , Table I) falls to the 
mght of the appropriate curve in figure 6. However, since the heating rate associated 
with this gas at this flight condition correlates well (as will be shown later) with the 
results for the other gases, it is possible that all of the correletions may be valid at 
higher flight speeds than indicated as figure 6. 

In figures 7 and 8 are shown a typical set of numerical solutions that have been 

obtained on an IBM 7090 digital computer for the boundary layer equations for the 

dissociating diatomic gases O 2> N 2 and H 2 , for a particular set of flight conditions. 

Figure 7 indicates the variation of gas temperature and tangential gas velocity across 

the stagnation point boundary layer for a vehicle travelling at 24, 000 ft. /sec. , at an 

2 

altitude such that the ambient pressure and temperature are 0.472 LB/FT. and 449°R 
respectively. The instantaneous surface temperature was taken to be 800°R for all 
three cases. Figure 8 shows the corresponding variation of gas composition through 
the layer as a function of r,, the non-dimensional distance measured from the surface. 

The resulting boundary layer distributions of temperature, tangential velocity 
and gas composition, in air and carbon dioxide, corresponding to the same gas 
conditions utilized for the pure diatomic gases, are shown in figures 9 to 11. By way 
of comparison, results for oxygen and nitrogen are also shown in figure 9. Here, the 
surface temperature was arbitrarily taken as 3500°R. 
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It is interesting to observe that whereas figures 8, 10 and 11 present the 


variation of equilibrium gas composition with distance from the surface in a layer of 
constant pressure, figures 3, 4 and 5 show the variation of equilibrium gas composition 
with gas temperature at constant pressure. Since it is seen in figures 5 and 7 that the 
gas temperature does not vary linearly with distance from the vehicle surface, clearly 
figures 8, 10 and 11 are essentially distorted versions of a set of curves which would 
look like figures 3, 4 and 5, if the former were replotted versus gas temperature. 

The boundary layer temperature and velocity profiles for a gas mixture 
composed of 50% nitrogen and 50% carbon dioxide by mass are shown in figure 12. 

The limiting cases of pure nitrogen and pure carbon dioxide are included for com- 
parison. The variation in gas composition of the nitrogen-carbon dioxide mixture, as 
a function of the non-dimensional distance within the boundary layer, is shown in 
figure 13. 

Note (see figure 12) that although the velocity profiles for these three cases 
are similiar, the temperature profiles are quite different. The nitrogen boundary layer 
is seen here to be characterized by a high temperature gradient at the surface. This is 
due to the fact that atomic nitrogen recombines at relatively high temperatures within 
the gas phase and this gas temperature is sustained rather close to the surface. In 
contrast with this, we note that in the carbon dioxide boundary layer, the chemical 
reactions (see figure 11) release energy further out in the gas phase, lowering the gas 
temperature and resulting temperature gradient at the surface. However, as will be 
shown later, the heat transferred to the surface by a carbon dioxide gas is greater than 
that transferred by a nitrogen gas (all other conditions being similiar), indicating, that 
in carbon dioxide, an appreciable amount of energy is transfered to the surface by the 
diffusion of chemical species, in addition to the energy transfered by conduction. 
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b) Heat Transfer and Skin Friction Correlations 

The aerodynamic heat transfer rate to the surface of the vehicle, in the absence 
of mass transfer or surface slip is given by: 


Q = ( K -P- - £ p.V.h.\ 

w l By i i 1 1 I 


(49) 
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In terms of the transformed variables, equation (49) becomes, when thermal diffusion 
is neglected: 
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where the modified Newtonian velocity gradient is given by: 
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and all other quantities appearing in equation (50) are obtained during the solution of 
the boundary layer equations. 

Upon introducing the following definitions for the Nusselt number and Reynolds 
number, based on surface conditions, 
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one obtains: 
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where: 
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( 55 ) 
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Two attempts have been made to correlate the effect of molecular weight of 
the various gas systems upon the dimensionless group Nu/vfte w , and these are shown 
in figures 14 and 15. When the numerical results from equation (50) are inserted into 
equation (54) for each of the eight gas systems treated here and these results are 
plotted versus flight velocity, one obtains the results shown in figure 14. As is shown, 
this method of correlation might be useful when only a single gas is involved, but when 
the results of other gases are included, no molecular weight correlation appears. When 
one plots the dimensionless heat transfer parameter Nu/\fRe as a logarithmic function 
of the normalized product of density and viscosity across the boundary layer, one 
obtains the results shown in figure 15. It is most interesting that a plot of Nu/vRe w 
versus o U / D U collapses the extensive numerical computations since here the 
effects of molecular weight, although certainly present, are not at all pronounced. It 
is the opinion of the authors that the data presented in figures 14 and 15 is not necessarily 
in the most useful form and that better methods exist for the correlation of hypersonic 
heat transfer. These will be shown below. 

Earlier work (e. g. , Ref 26) on aerodynamic heat transfer at hypersonic 
speeds suggests that another useful way of correlating the results would be with 
stagnation enthalpy and stagnation pressure. That is, one finds that the normalized 


aerodynamic heat transfer parameter Q ^ B is usually a linear function of the 

enthalpy difference (H - h w ) across the boundary layer, and when one plots the 
numerical results in this manner for the eight gaseous systems ti*eated here, one 
obtains the results shown in figure 16. Here, the molecular weight effect is most 


pronounced and is, in fact, separable. 

Thus, when one divides the aerodynamic heat transfer parameter Q, 


V >P~ 


by the enthalpy difference (H - h ), and plots this result versus the molecular 

G W 
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weight of the ambient undisturbed gas#, , one obtains the most useful result shown in 
figure 17. Here, it is seen that a correlation of all the numerical solutions to a high 
order set of coupled non-linear differential equations, governing the aerodynamic neat 
transfer in eight different chemically reacting gases, reduces to a single straight line. 

This straight line can be represented by the equation; (to within ±8%): 

, 1/2 


Q (R r /PJ 

w B e 
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e w 


(12. 0 + 0. 866 7ft J x 10 
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where the symbols have the following meaning and physical dimensions. 

BTU 
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heat transfer rate, 


ft. sec. 


(56) 
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H - 

e 

h - 

w 


vehicle nose rad.us, ft. 

stagnation oressure, atm. 

molecular weight of ambient gas 

stagnation enthalpy, BTU/lb. 

static enthalpy of gas mixture evaluated it 
surface temperature, BTU/lb. 


The results of other theoretical and experimental studies pertaining to 
hypersonic stagnation point heat transfer in several different gas mixtures are also 
included in figure 17. Item 6 on figure 17 indicates the range of the correlated theo- 
retical heat transfer parameter in air and is explained in greater detail in table IV. 
Hoshizaki's theoretical results for C0 2 are seen to fall slightly below the theoretical 
prediction of this present study. However, it is noted that there is favorable agree- 
ment between the various theoretical studies and experimental results. 
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Returning to equation (1), we are now in a position to relate the coefficient C 
to the first term on the right hand side of equation (56). We introduce the two hyper- 
sonic approximations 
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into equation (1) 
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and we obtain 
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where appropriate physical dimensions must be utilized in equation (59). Thus, for 
example, if one writes: 


< 6 ' 

Jd 

where all of the physical quantities in equation (60) have been defined immediately 
following equation (56), then one has: 


C = (6. 0 + 0.433 x 10" , 


lb. 

T 3/ 2 7 172 

it. sec. atm. 


On the other hand, if one writes: 



where the physical quantities have dimensions given by: 

P ro - ambient density, lb. /ft. 3 
Rg - nose radius, ft. 

V - flight speed, ft. /sec. 

00 

then the coefficient C has the value: 

C = (9. 18 + 0. 663 x 10 -10 , — TU sec ; 2 

ft. 3 lb 1/2 


(61) 
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Note that although equation (58) as written, applies specifically at the axi 
symmetric stagnation point of a body in a hypersonic stream, the result is actually 


considerably more general since one anticipates that 

6 


Q wV^B 

f (A) 

P e l/2 f(A) 

U%J 


-3 


= (12 + 0.866 jfcjx 10 (H e -h w ) 


(63) 


which applies over both two-dimensional (cylindrical) and axi-symmetric shapes 


including the effects of yaw, where the additional terms are defined as follows: 

6 is the Kroncker delta such that 


! 6 = 1: (axi-symmetric shapes) (64) 

6 = 0: (two-dimensional shapes) 


P is the local static pressure at the outer edge of the boundary layer expressed 
e 

in atmospneres, H e is the stagnation enthalpy in BTU/lb. , and h^ is the local wall gas 
enthalpy in BTU/lb. , and f (A) represents the effects of yaw and may be written: 

f (A) -- cos (65) 

to a good first approximation if the surface temperature remains relatively cool, 

(ref. 30). 

For a modified Newtonian pressure distribution 
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( 66 ) 


s 

and hence over the forward portion of the body where the first term on the right 
dominates, the heat transfer rate will decrease as one proceeds away from the forward 
stagnation point at a rate directly proportional to cos <p . 

It is noted, in passing, that for the same ambient conditions of density (p^) 
and flight velocity (V ro ) the heat transfer rate in pure carbon dioxide is 36% larger 
than in air. However, in a test facility, in which one may have greater control over 
To,,, Poo and V*, say, or in flight for which T^, Poo, and Voo are fixed, then the heat 
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transfer in a pure carbon dioxide atmosphere is 68% larger than in air, since for this 
situation the molecular weight effect shows up twice. 

Upon introducing the definition of the Stanton number and the shock Reynolds 


number 


H 



(67) 
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one obtains: 
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The skin friction coefficient may be defined as 
T 

C f = 2 
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where the aerodynamic shear stress at the surface of the vehicle is given as 



Hence, one also obtains: 
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Upon calculating the quantities C TT V R„ and P ) ( 5_) for all of the gaseous 

V c f4W 

systems and plotting the results versus the quantity V^/R^^-Jr) , one obtains a set 

of straight line correlations (see figure 18). It appears that one can therefore write: 
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( 75 ) 


It is noted that Reynolds' analogy states that 
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substituting equation (75) into equation (72) and using the relationship for the velocity 


given in equation (55) one obtains 



Clearly then, the results of figure 18 would appear to indicate that Reynolds' analogy 
applies to the same degree of accuracy as the Prandtl number of the various gases 
differs from unity. In figure 19 it is seen that the theoretically determined Prandtl 
number for a wide range of gases, lies in the range 0. 62 <: Pr <: 0. 74. 


CONCLUSIONS 

In this study of the aerodynamic laminar stagnation region heat transfer in 
eight different dissociating gas mixtures, a correlation of all of the numerical results 
to a high order set of coupled non-linear differential equations yielded a single straight 
line representation. That is, it was found that the aerodynamic heat transfer parameter 

"‘ s directly proportional to the enthalpy difference (H e ~h^) for each of the eight 
gases treated. Moreover, when the effects of molecular weight were correlated, ii 
was found that the aerodynamic stagnation point heat transfer parameter correlates 
linearly with the molecular weight of the undisturbed gas in the following manner. 
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It was also found possible to correlate the skin friction coefficient C^ for the 
eight gases treated, and it was shown that Reynolds' analogy provides a reasonable 
representation of the ratio of heat transfer to skin friction. 
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The results obtained here apply over a wide range of hypersonic flight speeds. 
Since more energy is required in order to dissociate C0 9 and H than air, the results 
apply directly at considerably higher flight speeds in a pure hydrogen atmosphere and 
at somewhat higher flight speeds in a pure carbon dioxide atmosphere than in Earth's 
atmosphere. It is expected that the results may also apply to the complex aerothermo- 
chemical environment produced by the gases flowing out of a rocket motor over a test 
model. 
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APPENDIX 

The transport properties can be described by the distribution function 
fj/ ^( r » t) which is the solution of the well known Maxwell- Boltzmann integro- 
differential equation. 



The right hand side of eq. A1 accounts for the collision processes occurring between 
the components of the system. The Maxwell-Boltzmann equation may be solved by 
means of the Enskog perturbation technique, which ultimately reduces the equations to 
linear integrals and finally to a set of linear algebraic equations. Chapman and 
Cowling (Ref. 31) and Hirschfelder, Curtiss and Bird (Ref. 23) have shown that the 
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coefficients of these algebraic equations may be expressed in terms of the definite 
integrals, 0(f ,s) . Recently, Cohen, Spitzer and Routly (Ref. 32), derived a modified 
equation which they considered more appropriate for inverse square forces, and 
treated the distant encounters by means of the Fokker-Planck equation. 

Following the Hirschfelder, Curtiss and Bird technique (Ref. 23) the omega 
integrals, from which the various transport coefficients may be computed, are ex- 


pressed as: 
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where the angle of deflection \ is a function of the relative velocity (g) and the impact 

parameter (b) expressed as: 
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where r is the lowest root for which the radical in the integrand vanishes. The 
m 

cross section (Q) is 
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The transport properties of the pure species may be calculated from the 
following equations. 
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which are the rigorous kinetic theory formulae for the viscosity and the binary diffusion 
coefficients respectively. Note that the self diffusion coefficients are also obtained from 
eq. A7 by a proper interpretation of symbols (i. e. j si). In order to calculate these 
properties one requires a knowledge of the collision diameter q, and the reduced collision 
integral O*. In the above equations M_ is the reduced mass of the binary interaction and 
o\j is the simple average of the collision diameters for the pure species i and j. 

Having a knowledge of the transport properties of the pure species and the 
composition of a gas mixture, one is then in a position to consider the transport 
properties of the gas mixture. 

The coefficient of viscosity for an n component gas mixture may be written 
in terms of the mole fractions, the viscosities of the pure species and the binary 
diffusion coefficients as follows: (Ref. 23) 
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Although the binary diffusion coefficients are symmetric (i. e. , 

multicomponent diffusion coefficients are not. This is apparent from the following 
definitions. (Ref. 23): 
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(All) 
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The frozen thermal conductivity of the mixture with the Eucken correction is given by 

(Ref. 33) K Eucken -K. 

Eucken = K . + l X. ^ 
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The elements of the above matrices are 
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The thermal conductivities of the pure species are given in terms of viscosities and 

specific heats as follows: 
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The quantities A*, and B* are defined in terms of the reduced omega integrals. 
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One, therefore, observes that in an n component gas there are n values of p,., 
n values of A. and — g— values of><9... This leads to one viscosity coefficient (^ Mix )» 
one thermal conductivity coefficient (K^^ ken ) and n 2 - n multicomponent diffusion 
coefficients (D.^). Finally, one is able to define the properties of the fluid mixture in 
terms of a single Prandtl number (Pr=M Mix C p /K^ ken ) and n 2 - n multicomponent 

Lewis numbers (L^-pC p D_/K^. x ). Note that we have utilized the following notation 
in the main text: 

„_„Eucken , , — - 

K “ K Mix. ; M= ^Mix. and C p = Cp Mix = ^ CjCp . 

i 

The variation of the Prandtl number with temperature within the boundary layer 
for the equilibrium dissociated mixtures considered in this study is shown in figure 19. 
It should be noted that the boundary layer pressure is different for each of the gases 
shown as indicated in table I. Note also that, since preferential diffusion effects are 
not neglected in this study, the relative concentrations of the chemical species for a 


given gas mixture are not alwhys in stochiometric proportions within the viscous layers. 
The magnitude of this effect may be assessed by comparison of the Prandtl numbers for 
air, nitrogen and oxygen given in figure 19. 
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NORMAL SHOCK AND STAGNATION POINT COMPARISONS 
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EQUILIBRIUM CONSTANTS, log _ K = a - b/T(°R) 
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TABLE IV 


COMPARISON OF CORRELATED THEORETICAL HEAT 
TRANSFER PARAMETERS IN AIR 

(calculations based on data appearing in ref, 5 below) 


INVESTIGATORS 

Reference 

<W P e> 1/2 . LB 

(H - H ) FT 3//2 SEC ATM 1//2 

G W 

SIBULKIN 

(1) 

0.034 - 0.042 

LEES 

(2) 

0.033 - 0.037 

ROMIG 

(3) 

0. 036 - 0.046 

FAY& RIDDELL (Le = 1.0) 

(4) 

0. 040 - 0. 043 

FAY & RIDDELL (Le = 1 . 4) 

(4) 

0. 045 - 0. 047 

SC A LA & BAULKNIGHT 

(5) 

0.035 - 0.043 
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igure 4. Equilibrium Dissociation of Carbon Dioxide 






Figure 6. Ambient Densities and Flight Velocities at Approximately One Per Cent Ionization 


Vqj = 24,000 FT/SEC.; P w = 0.472 LB/FT'S T C0 = 449°R 
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Figure 7. Temperature and Velocity Profiles (Diatomic Gases) 
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Variation of Atomic Species Through the Boundary Layer (Diatomic Gases) 
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Figure 12. Temperature and Velocity Profiles 
T co = 449°R, P = 0.472 lb/ft 2 
V ro = 24, 000 ft? sec, = 3500°R 
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Figure 15, Variation of the Nusselt Number with Product of Density and Viscosity 
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Figure 16. Variation of Correlated Heating Rates with Enthalpy Difference 
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Figure 17. Heat Transfer Correlation 






igure 19. Prandtl Numbers for the Equilibrium Dissociated Mixtu 
(P = 0.4715 lb/ft 2 , T = 449°R, V = 24, 000 ft/sec) 
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